Hidden parameters in open-system evolution unveiled by geometric phase 
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We find a class of open-system models in which individual quantum trajectories may depend on 
parameters that are undetermined by the full open-system evolution. This dependence is imprinted 
in the geometric phase associated with such trajectories and persists after averaging. Our findings 
indicate a potential source of ambiguity in the quantum trajectory approach to open quantum 
systems. 
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I. INTRODUCTION 

Closed quantum systems evolving deterministically un- 
der some Hermitian Hamiltonian is an idealized descrip- 
tion that at best approximates real laboratory exper- 
iments. In fact, all quantum systems undergo open- 
system effects induced by entanglement with environ- 
mental degrees of freedom; effects that may be detrimen- 
tal in various quantum information protocols in which 
coherence is an essential ingredient [l|. This feature has 
led to a revived interest in the theory of open quantum 
systems and how to deal with open-system effects by 
different types of error control to achieve error resilient 
quantum information processing [2j. 

Geometric and holonomic quantum computation, first 
conceived in Ref. Q and experimentally demonstrated 
in Ref. |4], is an approach to error control that has at- 
tracted considerable attention recently. In its simplest 
variant, it makes use of the Abelian geometric phase |5| 
to construct quantum logical gates acting on one or two 
quantum mechanical bits (qubits) [6j. These gates may 
be used to build quantum Boolean networks and may be 
combined with other error resilient methods to perform 
robust quantum computation [3, H|. The need to un- 
derstand the error resilience of geometric and holonomic 
quantum computation has led to proposals for the geo- 
metric phase of open quantum systems [9l-[l2| . 

Here, we examine the idea in Ref. [9J (see also Ref. 
|13j ) to associate geometric phases of individual quantum 
trajectories in quantum jump unravelings of Lindblad- 
type open-system evolution |14| . This approach involves 
only pure state geometric phases, which may relate to 
the geometry of the full open-system evolution by some 
averaging over trajectories [15J . 

The trajectory-based geometric phase simplifies the 
analysis of the robustness of geometric and holonomic 
quantum computation |r!J18| . The idea is that for weak 
influence of the environment, it suffices to consider the 
lowest order, no-jump trajectory to evaluate error re- 
silience. Here, we show that this geometric phase may 
in certain cases lead to different predictions regarding 
the resilience of geometric and holonomic quantum com- 
putation to open-system effects. This result indicates a 
potential source of ambiguity in the no-jump approach 



to analyze weak open-system effects. 

The problem of how to define open-system geometric 
phases by averaging over quantum trajectories has been 
addressed in Refs. [19ll20| . These works employ quantum 
state diffusion (QSD) [2f(, which is a form of stochastic 
unravelings of the Lindblad evolution consisting of con- 
tinuous, Brownian-type trajectories in state space. 

In Ref. [Hi, the averaged geometric phase associated 
with the full nonlinear form of the QSD equation [21[ 
was examined. It was found that this phase is not in- 
variant under unitary rotations L m —± ^ V mn L n of the 
Lindblad operators L m . On the other hand, Ref. [2fJ 
demonstrated that this noninvariance would disappear if 
the averaged geometric phase is instead associated with 
the linearized version of QSD [22|, provided the system 
starts in a pure state. Based on this result, it was claimed 
in Ref. [20| that the linearized QSD approach provides a 
uniquely defined geometric phase of open systems. Here, 
we demonstrate the existence of a class of Markovian 
open-system evolutions for which the linearized QSD ge- 
ometric phase may change under symmetry transforma- 
tions of the full Lindblad evolution. 

The outline of the paper is as follows. In the next sec- 
tion, we find symmetry transformations of a certain class 
of Markovian open-system evolutions which are not sym- 
metries of the corresponding no-jump trajectories. These 
transformations are shifts of the Lindblad operators, i.e., 
of the form L m —¥ L m — f m (t)l . Here, f m (t) are arbitrary 
complex- valued functions and are hidden parameters in 
the sense that they do not affect this class of Markovian 
evolution models. In Sec. IIIIl this result is illustrated by 
an explicit calculation of the no-jump geometric phase 
for a dephasing qubit (spin— i) being exposed to a static 
magnetic field. The geometric phase for stochastic un- 
ravelings in the form of the linearized QSD equation is 
analyzed in Sec. IIVI The paper ends with the conclu- 



II. SHIFT SYMMETRIES OF OPEN-SYSTEM 
EVOLUTIONS 

We consider Markovian evolution of open quantum sys- 
tems governed by the Lindblad equation (H = 1 from now 



on) [22 



p{t) = -t[H(t),p(t)}+\J2[L m p(t)Ll 
= -i[H{t),p{t)]+\£{p(t)). 



(1) 



Here, L m are dimensionless Lindblad operators that 
model the influence of the environment on the system 
evolution. For simplicity, we shall assume that L m are 
time-independent. The parameter A > controls the 
strength of the open-system effect, such that A = cor- 
responds to unitary closed system evolution. 

The Lindblad equation obeys certain symmetries; an 
apparent one is the independence of choice of zero point 
energy corresponding to the transformation H(t) — > 
H(t)—h(t)l, where h{t) is real valued and 1 is the identity 
operator. Another general type of symmetry corresponds 
to the transformation L m — > V^ V mn L n , where V mn is 
an arbitrary unitary matrix [20j . One may check that this 
transformation leaves the Lindblad equation unchanged 
and thus will not affect the state p{t) of the system. 

The quantum jump unraveling is defined by dividing 
the evolution given by Eq. ([1]) into small time steps At. 
In the At — > limit, this procedure leads to quantum tra- 
jectories in state space consisting of smooth deterministic 
parts interrupted by random jumps, generated by jump 
operators proportional to L m . For a pure initial state ipo, 
these trajectories reside in projective Hilbert space V{K) 
formed by rays of the system's Hilbert space H. These 
rays are equivalence classes consisting of vectors that dif- 
fer by multiplication of nonzero complex numbers. As 
shown in Ref. [9j], one may associate a pure state ge- 
ometric phase to each such trajectory. Here, we focus 
on the geometric phase of no-jump trajectories. Such a 
trajectory is the projection onto V{T-L) of the continuous 
Hilbert space path 

[0,T] 3t^\Ht)) =Te-^o"^^ dt '\^o) (2) 

with T time ordering and 

H(t)=H(t)- l -\Y,LlL m (3) 



a non-Hermitian effective no-jump Hamiltonian. The 
corresponding no-jump geometric phase acquired on the 
time interval [0, T] reads [9| 



lnj = arg(^(0)|V(T)) 



(m\H(t)m)) 
m)m)) 



dt. (4) 



Note that "f n j is a property of a path in "P("H) as it is 
invariant under the transformation \ip(t)) — > c(t)\i]j(t)) 

together with H(t) — > H(t) + i^ In t^t, where c(t) is a 

nonzero complex number for all t £ [0, T]. 



It is straightforward to check that the no-jump path 
\ip{t)} and the corresponding geometric phase j n j are un- 
affected by the above-mentioned unitary rotation L m —t 
X^n VmnLn- But there may be other symmetries that ap- 
ply only to certain classes of open systems. We focus on 
the symmetry related to the shifts L m — > L m — f m (t)l 
where f m (t) is in general complex- valued. Such shifts 
induce the transformations 



H(t) -► K(t) = H(t) 
C -> C. 



^Y.^*r, l {t)L m ~f m (t)Ll) 1 

m 

(5) 



Thus, they result solely in an extra term in the Hamilto- 
nian part of Eq. ([TJ . This implies that the Lindblad evo- 
lution is unchanged under the shifts of L m if all f* n (t)L m 
are Hermitian. In such a case, f m (t) are said to be hid- 
den parameters of the full open-system evolution. On the 
other hand, the no-jump Hamiltonian transforms as 



H(t)^K(t) = H{t)+ l -\Y J {f m {t)Ll 

m 
+ f m (t)L m -\frn(t)\ 2 l) (6) 



with H(t) the no-jump Hamiltonian in Eq. (j3|). The 
transformed no-jump Hamiltonian K(t) may be nontriv- 
ially different from H(t) even for Hermitian /^(i)L m . 
Thus, for shifts L m — > L m — f m (t)l such that /^(i)L m 
are Hermitian, the Lindblad evolution is unchanged but 
the deterministic no-jump evolution may undergo a non- 
trivial change originating from the anti- Hermitian contri- 
butions |A (/ m (t)L+, + f^(t)L m ) to K(t). In this case, 
the no-jump trajectories may depend on the parameters 
f m , which are hidden in the full open-system evolution. 
We note that this result applies to any smooth portion 
of a quantum trajectory, i.e., trajectories that contain 
one or several jumps share with the no-jump trajectories 
the same kind of behavior under shifts of the Lindblad 
operators. 

The f m dependence may be interpreted as an manifes- 
tation of a continuous monitoring of the environment in 
the presence of a specific form of system-environment in- 
teraction. To see this explicitly, let us consider a unitary 
representation model for the system-environment evolu- 
tion during the time interval [t,t + St], where St is the 
finite time resolution for measuring projectively the en- 
vironment in some orthogonal basis {|0 e ), |m e )}. We as- 
sume that St and A are much smaller than the typical 
energy shift associated with H(t). Under this assump- 
tion, the change in the system-environment state can be 
described by a unitary map U(t, t + St; {f m (t)}) with the 



effect 



|O e )|V(t)) -> U(t,t + 5t;{f m (t)})\O e )m)) 

= |o e ) (i - iK(t)6t) [V(*)> 

+VA^J] |m e ) (L m - / m (t)i) Mt)}. 

m 

(7) 

Here, we have assumed that the environment is prepared 
in the pure state |0 e ) and we have taken 

(0 e \U(t,t + St;{f m (t)})\0 e ) = l-iK(t)6t, 
(m e \U{t,t + St;{f m (t)})\O e ) = Vm(L m 

-fm{t)V) (8) 

to the first order in St and vXSt. Thus, the shifts 



L m — > L m — f m (t)l would correspond to engineer- 
ing the system-environment interaction so that Eq. 
([8]) is satisfied. Evidently, the jump operators are 
>/\5i (L m — f m (i)i) . The no-jump trajectory [0,T] 9 

£ -> |-0(i)) = Te'^o ^(*') d *'|V'o} is realized with proba- 
bility (V'(T)I^(T)) by verifying that no change has oc- 
curred in the environment, and repeating up to time T . 

The operators Fo(t) = 1 — iK(t)5t and F m {t) = 
\/XSt \L m — f m (t)l] in Eq. ((SJ) constitute a set of Kraus 
operators that represent a completely positive map of 
system states from t to t + St. Provided all f m (i)L m are 
Hcrmitian, there is a unitary matrix W(t) that relates 
this Kraus representation with the original one consisting 
of E (t) = 1 - iH{t)8t and E m (t) = V\5tL m . Explicitly, 
we may write F M (i) = J2 V W tiy (t)E v (t), /x, v = 0, 1, . . ., 
with 



W(t) 



A-|A«E m l/m(*)| 2 

-Vxsifi(t) 

-VXSif 2 (t) 

-Vxsify(t) 
V 




(9) 



which can be checked to be unitary to lowest order in St 
and yXSt. The existence of such a unitary matrix demon- 
strates (l| that the two maps are physically identical if 
all f^{t)L m are Hermitian; a result that is consistent 
with the invariance of the Lindblad equation under the 
shifts L m — > L m — f m (t)l for this type of open-system 
evolution. 




FIG. 1: (Color online) The no-jump phase 7 nj as a function 
of the open-system strength A for the shift values / = 0, 0.2, 2. 
The initial state is taken to be pure on the equator (9o = 5 ) 
of the Bloch sphere. The horizontal curve for / = demon- 
strates the A independence found in Refs. [y, [lj. For / ^ 0, 
7„j depends strongly on A. 



III. GEOMETRIC PHASE 

We next show that the previous result may have conse- 
quences for the geometric phase of a no-jump trajectory. 
We find an explicit physical example where a shift pa- 
rameter / is imprinted in the no-jump geometric phase, 
although the corresponding open-system evolution is /- 
independent. 

Consider a qubit (spin— h) prepared in the pure state 

IV'o) = cos (i^o) |0) + sin (5^0) |1), exposed to a static 
magnetic field in the z direction and to dephasing of 
strength A. This may be modeled by a Hamiltonian 
H = %<t z and a single Lindblad operator L = er 2 , 
with w the precession frequency and a z = |0)(0| — |1)(1| 
the z component of the standard Pauli operators. By 
using Eqs. ([5]) and ©, we find the transformations 
H -> K = H - AIm(/)cr z , C -> C, and H -> K = 
H + i\Re{f)a z ~ |A|/| 2 1, under the shift L -> L - /I, 
where / is assumed to be time- independent for simplicity. 
Thus, the Lindblad evolution is unaffected by the shift if 
/ is real-valued. On the other hand, the geometric phase 
of the no-jump trajectory 

[0,T]3t^\m)=e- lRt m (10) 

for a quasi-cyclic path over the time interval [0, 2tt/uj] 



with real- valued /, takes the form 



Inj 



n + If\ ln { e ~ 
fe -4?A s in 2 ^ 



'-fX 2"0 

J COS — - 



(11) 



which is explicitly /-dependent. 

The geometrical reason for this /-dependence can be 
seen by looking at the Bloch sphere polar angle 0, 
which becomes time-dependent if / 7^ 0. Explicitly, 
by evaluating the right-hand side of Eq. (fTUf we ob- 
tain tan[0(i)/2] = e~ 2 ^ A 'tan(# /2) and azimuthal angle 
tp(t) = (fo + tot, which correspond to a spiralling motion 
toward the north (south) pole of the Bloch sphere for 
/ > (/ < 0) and all O ^ tt (0 o # 0). Furthermore, 
one may check that j n j in Eq. (jlll) converges to the ex- 
pected — 7r(l — cos 6*o ) (minus half the solid angle enclosed 
on the Bloch sphere) in the /A — > limit. The nontrivial 
/ dependence is illustrated in Fig. [TJ The resilience to 
dephasing of the geometric phase of the no-jump trajec- 
tory found in Refs. [9|, [l7[ corresponds to the case / = 0. 
However, as our calculation shows, any nonzero / would 
predict -f n j to be A-dependent and thus be affected by this 
kind of open-system effect. For small /A, this dependence 
is linear, which may be seen by expanding the geomet- 
ric phase around the closed system expression, leading to 
the lowest order correction 2tt 2 L- sin 2 a . 

We may show that the no-jump evolution of dephasing 
for / 7^ is equivalent to decay of the precessing qubit. 
Consider the evolution generated by the Lindblad oper- 
ator L- = a x — i<Jy, which corresponds to decay toward 
the south pole of the Bloch sphere with some strength A', 
say. The no-jump curve is determined by the effective no- 
jump Hamiltonian H' — ^a z —iX'a z —iX'l, where we have 
assumed the Hamiltonian H = ^a z . If A' = — /A, then 

the no-jump Hamiltonian H' generates the same curve 
in P(7i) as H in the dephasing model (the correspond- 
ing Hilbert space curves differ only by multiplication of 
a nonzero complex number). 

It is instructive to compare the preceding dephasing 
example with the mixed state geometric phase j[P] pro- 
posed in Ref. [l(| for real- valued /. Since j[V] is based 
directly on the kinematics p(t) of open-system evolution, 
it follows immediately that "/[P] is /-independent and 
therefore experimentally testable |24| . It was further- 
more found that j[P] for precession around the z axis 
is resilient to dephasing only if the initial state lies on 
the equator of the Bloch sphere. This particular form of 
resilience has been reported in a recent experiment [251 ] 
with polarized ultracold neutrons exposed to dephasing 
noise. 

For one-qubit systems, f m (t)L m can be chosen to be 
Hcrmitian for all open-system models that include Lind- 
blad operators that are linear combinations of the Pauli 
operators with real coefficients, such as dephasing and 
depolarization. On the other hand, it should be stressed 
that if no f m (t) exist such that f^(t)L m becomes Hermi- 



tian, then the shift in fact corresponds to a new Hamil- 
tonian that may cause a different evolution p(t). One 
such example is spontaneous decay of a qubit, for which 
one cannot find a nonzero / such that /*£- becomes 
Hcrmitian. Indeed, the shift L_ — » L_ — /l induces the 
Zeeman term A \hn( f)<j x + Re(/)(jj,] in the Hamiltonian, 
in the case of spontaneous decay. Another example is 
a spin— i system interacting with a quantized light field 
and subjected to a linear loss of photons. This loss may 
be modeled by the non-Hermitian photon annihilation 
operator L — a. In this case, the shift L — > L — /l yields 
the extra term \/2 [— Im(/)x + Re(/)p] in the Hamilto- 
nian, where x — (a + a T )/V2 and p = (a — a')/(iv2) are 
the 'position' and 'momentum' operators, respectively, of 
the quantized light field. It is reasonable to expect that 
any robustness in the geometric phase found in these two 
latter types of systems would signal an error resistance 
also at the level of the full Lindblad evolution. Appli- 
cations of the quantum jump unraveling to these model 
systems have indeed shown a nontrivial error resilience 
[la . [26f , results that support the usefulness of the geo- 
metric phase for robust quantum computation. 



IV. STOCHASTIC UNRAVELINGS 

Stochastic unravelings in the form of quantum state 
diffusion (QSD) [21| consist of continuous, Brownian-like 
quantum trajectories whose average coincides with the 
full Lindblad evolution. The geometric phase of such tra- 
jectories arising from nonlinear 21 [ and linear [22j ver- 
sions of the QSD equation has been considered in Refs. 
19] and 20], respectively. Reference [19] considered the 



phase transformation L ri 



e %Xm L m and found a non- 



trivial Xm dependence in the geometric phase for the non- 
linear QSD evolution. In Ref. [2(|, it was demonstrated 
that the averaged geometric phase a g associated with the 
linearized evolution is invariant under unitary rotations 
L m -^ J2n V-mnLm provided the system starts in a pure 
state. Here, we examine the behavior of this geometric 
phase under the shifts L m — > L m — f m (t)l and show that 
a g may depend on f m , also when f m is hidden in the full 
open-system evolution. 

The linearized QSD equation reads 



|ety) 



1 



-iH(t)dt--Xj2 L l L m 

m 

VxY.LmdwJ^), (12) 



where w m are complex Wiener processes with respect to a 
probability measure Q. There is a mean Eq over Q such 
that E,Q[dw m ) = Eq[dw m dw m '} = 0; EqldWmdw^,] = 
5mm 1 dt. This guarantees the properly renormalized av- 
erage of any measurable quantity to coincide with the 
expectation value with respect to p(t). Following Ref. 
[20J . the averaged geometric phase a g with respect to 



the probability measure Q is taken to be 



argE Q [(0 o |<XT))] 



Tr [p(t)H(t)} dt. (13) 



The second term on the right-hand side of this expression 
depends only on the full state p(t) and would therefore 
be unaffected under all symmetry transformations of the 



Lindblad equation. To show the noninvariance of a g un- 
der shifts of the Lindblad operators, it is therefore suffi- 
cient to show that the first term may be f m dependent. 
We demonstrate this by an example, again the dephasing 
qubit model with real- valued and time-independent shift 
parameter / and Hamiltonian H = %o z . For initial state 

|0o) = cos (5^0) |0) + sin (56*0) |1)> we obtain 



J 



argE Q [((/>o|(£(T))] = arg(0 o | exp 



1 

— uj ■ 
2 



ifX To z 



|</i>o) = — arctan 



tanh(/AT) + cos 6 
1 + tanh(/AT) cos 6 



tan 



~2 



(14) 



r 



Thus, argEQ[(0o|0(T 1 ))] is / dependent if uT 7^ n7r, n 
integer, and cos#o ^ ±1. Thus, it follows that the aver- 
aged geometric phase a g associated with the linearized 
QSD evolution may depend on the hidden parameter /. 



V. CONCLUSIONS 

We have demonstrated the existence of Markovian 
open-system evolutions for which the associated no-jump 
quantum trajectories may depend on parameters that are 
undetermined by the full open-system evolution. We 
have found conditions for this situation to occur and 
have identified the origin of this dependence in terms 
of continuous monitoring of the system's environment. 
Furthermore, we have explicitly demonstrated how such 
a hidden parameter can be unveiled by the geometric 



phase of an individual quantum trajectory for a dephas- 
ing qubit. The realization of the geometric phase for 
single quantum trajectories requires explicit engineering 
of the system-environment interaction; a feature that is 
shared by the mixed state geometric pha ses for com- 
pletely positive maps proposed in Ref. [271 ] . Finally, we 
have demonstrated that the averaged geometric phase in- 
troduced in Ref. [20| of the linearized QSD model shows 
a similar dependence on hidden parameters. Thus, it 
remains open whether a well-defined open-system geo- 
metric phase based upon quantum trajectories exists. 
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